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Kolomeisky, Newman, Straley and Qi Reply: 

Bhaduri and Sen (BS) § criticize our theory of low- 
dimensional inhomogeneous Bose liquids B on two points 
which we refute below. 

1. BS object to the way our energy functional depends 
on the gradient of the density, and propose two alterna- 
tives. As they admit, the first of these (no dependence) 
fails to describe what happens at a hard wall. It is also 
inferior to our theory in representing what happens near 
a classical turning point (see Figure 1 of our paper ||), 
and is unable to describe the important quantum me- 
chanical effect of particle tunnelling through a classical 
barrier, which is present in our theory. 

Their second choice is even worse: the term — {p') 2 /p 
means that gradients lower the energy: the ground state 
is unstable against spontaneous waviness. The corre- 
sponding nonlinear Schrodinger equation docs not have 
physically sensible solutions. 

Figure 1 compares the behavior near x — 
of the exact density function for a dark soli- 
ton § n(x) = (2N + 1)/2L - sin[27r(2iV + 
l)x/L]/[2Lsbx(27rx/L)] with the approximate expression 
n apP rox(x) = (2N/L)[3cot}i 2 (kFx) — l]^ 1 which is Eq. 
(12) of our Letter, with (3 = (however, please note that 
fci? = 2ttN/L for this example, because we only fill every 
other single particle state). 



A = p(n 



x-l/2 



h(y/p(r)), e = 6 Q {r) + \^y\ (2) 




FIG. 1. Exact and approximate density profiles near a soli- 
ton. The exact solution has additional short wavelength struc- 
ture not captured by the long-wavelength theory. 

The alternate forms for the energy functional proposed 
by BS cannot describe this stationary soliton at all. The 
density gradient terms in our theory are appropriate and 
necessary. 

Gray soliton solutions of the system of free fermions 
have been generated numerically fl although their exact 
analytical form is unknown. Here our theory is the only 
source of analytical results. 

Another example which verifies the starting point of 
our theory and can be compared with the exact re- 
sults concerns self-similar solutions whose existence was 
pointed out in our Letter 0. Consider a dilute system 
of particles placed in a harmonic trap of dimensionless 
frequency w. In the dimensionless variables of our Letter 
it is described by 



2l3rf = -d 2 f + [|/| 4 + W V]/ 



The self-similar solution / = Ae l9 derived from (|l|) has 
the form 



where the functions p(r) and h(v) obey the equations 



fp 

dv 2 



-w p+ — 

P 6 



= -Sh + h + ryv h 



(3) 
(4) 



where 7 and 5 are arbitrary constants. Eq.(H) for the scal- 
ing function h(v) has localized solutions only for 6 > 
and 7 > 0: for 7 = an explicit analytical solution 
can be written down, while for 7 > 0, (Q) has the same 
functional form as the equation we encountered in deter- 
mining the density profile in the harmonic trap [cf. Eq. 
(6) of the Letter for V = muj 2 x 2 /2]. 

The dynamics of the length scale p(r) can be under- 
stood by viewing (|) as a fictitious classical mechanics 
problem in the potential U = w 2 p 2 /2 + -f/2p 2 . This anal- 
ogy implies that an initially localized cloud of particles 
in free space (w = 0) expands asymptotically in a ballis- 
tic fashion: p(r) ~ t. These conclusions are in complete 
agreement with the results of a full quantum- mechanical 
calculation Q of expanding cloud of free fermions. In 
the presence of the confining potential (w ^ 0) the scale 
p(r) oscillates between maximum and minimum values: 
for 7 = the dynamics of p is the same as for a harmonic 
oscillator of frequency w. 

2. BS state that current experimental systems are not 
dilute. However, appropriate experiments will be possi- 
ble in the future (their feasibility has been discussed ||) 
to which our results are relevant. Our approach can pro- 
vide a useful insight even where it is not quantitatively 
correct. 
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